subroutine oleic_acid_init()
use molecule
implicit none
   !-------------------------------------------------------------------------
   natom    = 54
   neighmax = 4
   molname  = "OleicAcid"
   !
   if ( allocated(neigh   )) deallocate( neigh   )
   if ( allocated(element )) deallocate( element )
   if ( allocated(molpos  )) deallocate( molpos  )
   allocate( element(natom), neigh(neighmax, natom) ) !, molpos(3,natom) )
   !
   element( 1) = "HO"
   element( 2) = "OH"
   element( 3) = "C "
   element( 4) = "O "
   element( 5) = "CT"
   element( 6) = "HT"
   element( 7) = "HT"
   element( 8) = "CT"
   element( 9) = "HT"
   element(10) = "HT"
   element(11) = "CT"
   element(12) = "HT"
   element(13) = "HT"
   element(14) = "CT"
   element(15) = "HT"
   element(16) = "HT"
   element(17) = "CT"
   element(18) = "HT"
   element(19) = "HT"
   element(20) = "CT"
   element(21) = "HT"
   element(22) = "HT"
   element(23) = "CT"
   element(24) = "HT"
   element(25) = "HT"
   element(26) = "CM"
   element(27) = "HC"
   element(28) = "CM"
   element(29) = "HC"
   element(30) = "CT"
   element(31) = "HT"
   element(32) = "HT"
   element(33) = "CT"
   element(34) = "HT"
   element(35) = "HT"
   element(36) = "CT"
   element(37) = "HT"
   element(38) = "HT"
   element(39) = "CT"
   element(40) = "HT"
   element(41) = "HT"
   element(42) = "CT"
   element(43) = "HT"
   element(44) = "HT"
   element(45) = "CT"
   element(46) = "HT"
   element(47) = "HT"
   element(48) = "CT"
   element(49) = "HT"
   element(50) = "HT"
   element(51) = "Ct"
   element(52) = "Ht"
   element(53) = "Ht"
   element(54) = "Ht"
   !
   neigh(:, 1) = (/ 2, 0, 0, 0 /)
   neigh(:, 2) = (/ 1, 3, 0, 0 /)
   neigh(:, 3) = (/ 2, 4, 5, 0 /)
   neigh(:, 4) = (/ 3, 0, 0, 0 /)
   neigh(:, 5) = (/ 3, 6, 7, 8 /)
   neigh(:, 6) = (/ 5, 0, 0, 0 /)
   neigh(:, 7) = (/ 5, 0, 0, 0 /)
   neigh(:, 8) = (/ 5, 9,10,11 /)
   neigh(:, 9) = (/ 8, 0, 0, 0 /)
   neigh(:,10) = (/ 8, 0, 0, 0 /)
   neigh(:,11) = (/ 8,12,13,14 /)
   neigh(:,12) = (/11, 0, 0, 0 /)
   neigh(:,13) = (/11, 0, 0, 0 /)
   neigh(:,14) = (/11,15,16,17 /)
   neigh(:,15) = (/14, 0, 0, 0 /)
   neigh(:,16) = (/14, 0, 0, 0 /)
   neigh(:,17) = (/14,18,19,20 /)
   neigh(:,18) = (/17, 0, 0, 0 /)
   neigh(:,19) = (/17, 0, 0, 0 /)
   neigh(:,20) = (/17,21,22,23 /)
   neigh(:,21) = (/20, 0, 0, 0 /)
   neigh(:,22) = (/20, 0, 0, 0 /)
   neigh(:,23) = (/20,24,25,26 /)
   neigh(:,24) = (/23, 0, 0, 0 /)
   neigh(:,25) = (/23, 0, 0, 0 /)
   neigh(:,26) = (/23,27,28, 0 /)
   neigh(:,27) = (/26, 0, 0, 0 /)
   neigh(:,28) = (/26,29,30, 0 /)
   neigh(:,29) = (/28, 0, 0, 0 /)
   neigh(:,30) = (/28,31,32,33 /)
   neigh(:,31) = (/30, 0, 0, 0 /)
   neigh(:,32) = (/30, 0, 0, 0 /)
   neigh(:,33) = (/30,34,35,36 /)
   neigh(:,34) = (/33, 0, 0, 0 /)
   neigh(:,35) = (/33, 0, 0, 0 /)
   neigh(:,36) = (/33,37,38,39 /)
   neigh(:,37) = (/36, 0, 0, 0 /)
   neigh(:,38) = (/36, 0, 0, 0 /)
   neigh(:,39) = (/36,40,41,42 /)
   neigh(:,40) = (/39, 0, 0, 0 /)
   neigh(:,41) = (/39, 0, 0, 0 /)
   neigh(:,42) = (/39,43,44,45 /)
   neigh(:,43) = (/42, 0, 0, 0 /)
   neigh(:,44) = (/42, 0, 0, 0 /)
   neigh(:,45) = (/42,46,47,48 /)
   neigh(:,46) = (/45, 0, 0, 0 /)
   neigh(:,47) = (/45, 0, 0, 0 /)
   neigh(:,48) = (/45,49,50,51 /)
   neigh(:,49) = (/48, 0, 0, 0 /)
   neigh(:,50) = (/48, 0, 0, 0 /)
   neigh(:,51) = (/48,52,53,54 /)
   neigh(:,52) = (/51, 0, 0, 0 /)
   neigh(:,53) = (/51, 0, 0, 0 /)
   neigh(:,54) = (/51, 0, 0, 0 /)
   !
   !-------------------------------------------------------------------------
end subroutine oleic_acid_init
